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TECHNICAL PUBLICATION 


USING THE INFLECTION POINTS AND RATES OF GROWTH AND DECAY 
TO PREDICT LEVELS OF SOLAR ACTIVITY 

1. INTRODUCTION 


Sunspot cycles are conventionally described using 12-mo moving averages (12-mmas) of 
monthly mean sunspot number, where (since 1981) the official number is now determined by the 
Sunspot Index Data Center <http://sidc.oma.be/index.php3>, 1 located at the Royal Observatory 
of Belgium (in Brussels). Formerly, it was maintained at the Swiss Federal Observatory in Zurich, 
Switzerland. The minimum value of the 12-mma of monthly mean sunspot number is called the 
sunspot minimum amplitude (Rm) and its occurrence denotes the epoch of Rm ( E(Rm )). Likewise, 
the maximum value of the 12-mma of monthly mean sunspot number is called sunspot maximum 
amplitude ( RM) and its occurrence denotes the epoch of RM ( E(RM )). In reality, sunspot mini- 
mum and maximum are better pictured as being broad intervals of time of several years in length 
when sunspot numbers are predominantly lower and higher, respectively, rather than being specific 
instances in time. While this is true, it is convenient to use the minimum and maximum values and 
their epochs of occurrence to describe the general characteristics of a sunspot cycle, such as its rela- 
tive size, ascent duration (ASQ, descent duration ( DES ), and cycle length (called period (PER)). 

Following Rm, the sunspot number gradually increases in value, typically reaching a maxi- 
mum rate of growth about 1-2 yr after E(Rm) and attaining RM about 3-5 yr after E(Rm). The point 
of maximum positive rate of growth in sunspot number represents the ascending inflection point, 
and its numerical value has proven to be useful for estimating the later occurring RM. 2 10 Likewise, 
during the declining portion of the sunspot cycle, there is a second inflection point (the descending 
inflection point), occurring typically about 6-7 yr after Rm, which seems to be related to the period 
of the ongoing sunspot cycle. 2 

The purpose of this Technical Publication is to reexamine the role of the inflection points and 
rates of growth and decay in sunspot cycle prediction on the basis of the behaviors of cycles 12-23, 
the most reliably determined sunspot cycles. 11-13 



2. RESULTS AND DISCUSSION 


2.1 An Overview of Cyclic Behavior (Cycles 12-23) 

Figure 1 displays the mean cycle curve, based on epoch analysis and cycles 12-23 (the smoothed 
thick line) for elapsed time in months (?) from E(Rm), t- 0-132 mo. The thin lines plotted above and 
below the mean cycle curve represent the upper and lower 90-percent prediction limits about each 
monthly value, respectively. The occurrence dates of E{RM) relative to E(Rm), which determines 
the ASC for each of the cycles, are noted across the top, and the actual RM values for the cycles are 
noted to the right. The respective PER for each of the cycles is noted near the bottom right, where 
cycle 23’s PER is not shown since it remains ongoing. It should be noted, however, that the National 
Oceanic and Atmospheric Administration (NOAA) Solar Cycle 24 Prediction Panel 14 has predicted 
cycle 24’s official onset to occur about March 2008 ± 6 mo, inferring that cycle 23’s PER will measure 
about 141 ±6 mo; hence, it will be a cycle of longer cycle length, like cycles 12-14 and 20 (the first 
new cycle spot was reported 1 in January 2008). 

Inspection of figure 1 suggests that cycles 12-23 might be loosely grouped either into three 
arbitrary groups based on relative size (large-amplitude cycles (18, 19, 21, and 22), average-ampli- 
tude cycles (15, 17, 20, and 23), and small-amplitude cycles (12, 13, 14, and 16)) or possibly into two 
groups based on cycle length (short-period cycles (15-19, 21, and 22) and long-period cycles (12-14, 
20, and 23)). 15 Table 1 gives averages and standard deviations (in parentheses) for Rm, RM, ASC, 
and PER for each of these groupings. It is found that large-amplitude cycles tend to rise more quickly 
to RM (shorter ASC), be of shorter PER, and have higher Rms than average-size cycles, which in 
turn have higher/shorter values as compared to small-amplitude cycles. Likewise, short-period cycles 
tend to rise more quickly to RM and to have higher Rms and RMs than long-period cycles. 

Figure 2 shows the mean cycle curve of the month-to-month rate of change in the smoothed 
monthly mean sunspot number (R), based on epoch analysis and cycles 12-23 (the smoothed thick 
line) for elapsed time in months from E(Rm), t- 0-132 mo. The occurrence dates for each cycle’s 
maximum positive rate of change during the ascending portion of the sunspot cycle are across the 
top, and the actual values of the maximum positive rate of change during the ascending portion of 
the sunspot cycle are to the upper right. The occurrence dates for each cycle’s maximum negative 
rate of change during the descending portion of the sunspot cycle are at the bottom, and the actual 
values of the maximum negative rate of change during the descending portion of the sunspot cycle 
are to the lower right. 

Based on the aforementioned groupings (amplitude and period), it is found that large-ampli- 
tude cycles tend to have higher rates of change both during the ascending (9. 5(1. 4)) and descending 
(-7.5(1 . 1)) portions of the sunspot cycle as compared to cycles of average size (6.7(1 .3) and -5. 2(0. 3)) 
and cycles of smaller amplitude (4.7(1. 1) and -4. 7(0. 9)). Also, cycles of shorter period tend to have 
higher rates of change during both the ascending (8. 5(2.0)) and descending (-6. 5(1. 5)) portions of 
the sunspot cycle as compared to cycles of longer period (4. 9(1.0) and -4.7(0. 8)). The timing of the 
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1. Variation of R for elapsed time in months t = 0-132 mo from E(Rm). 


Table 1 . Average parametric values for selected groupings of cycles. 


Group 

Cycles 

<Rm> 

<RM> 

<ASO 

<PER>* 

Amplitude: 


Large 

18, 19, 21,22 

8.9 (4.2) 

169.0 (22.1) 

40.5 (5.4) 

121.8(4.2) 

Average 

15, 17, 20,23 

5.6 (3.8) 

114.0 (7.3) 

46.8 (2.6) 

128.3 (10.4) 

Small 

12-14, 16 

3.9(1.7) 

76.2 (9.8) 

52.8 (6.4) 

134.3 (9.3) 

Period: 


Short 

15-19, 21,22 

6.6 (4.3) 

139.8(41.4) 

44.1 (7.1) 

121.9 (3.3) 

Long 

12-14, 20, 23 

5.5 (3.3) 

91.6 (23.8) 

50.2 (5.6) 

139.0 (2.9) 


* txcludes Cycle 2ti 


maximum positive rate of change in R (T/hnaxpos) and the maximum negative rate of change of 
R (Aifrnaxneg) are found to span 14-42 and 51-83 mo, respectively, and to average 27.8 (10.1) and 
71.4 (9.5) mo, respectively. 

Figure 3 depicts the 12-mma of the month-to-month rate of change shown in figure 2. The 
thick line is the mean curve and the thin lines above and below the mean curve represent the 90- 
percent prediction limits. Smoothing removes a considerable amount of the choppiness that appears 
in figure 2. The occurrences of the 12-mma of T/hnaxpos are across the top, the individual values of 
the 12-mma of T/hnaxpos are to the upper right, the occurrences of the 12-mma of ARi naxneg are 
across the bottom, and the individual values of the 12-mma of Tihnaxneg are to the lower right. 

Based on the aforementioned groupings (amplitude and period), it is found that large-ampli- 
tude cycles tend to have higher rates of change both during the ascending (6. 8(1. 3)) and descending 
(-4. 6(0. 6)) portions of the sunspot cycle as compared to cycles of average size (3. 9(0. 3) and -3.0(04)) 
and small-amplitude cycles (2. 7(0. 6) and -2. 3(0. 5)). Also, cycles of shorter period tend to have higher 
rates of change during the ascending (5. 5(1. 9)) and descending (-3. 8(1. 2)) portions of the sunspot 
cycle as compared to cycles of longer period (3. 1(0.8) and -2. 6(0. 7)). The timing of the 12-mma of 
ARmaxpos and T/hnaxneg is found to span 14-41 and 66-96 mo, respectively, and is found to aver- 
age 23.7 (7.8) and 76.8 (8.3) mo, respectively. 

Figure 4 depicts the cyclic variation of the previously discussed parameters. The horizontal 
line in each subpanel is the mean. The standard deviation ( sd ) is also given. It is noticeable that 
cycles of late (cycles 18-23) have higher averages of Rm, RM, ARmaxpos, ARmaxneg, 12-mma of 
ARmaxpos and 12-mma of Tifrnaxneg as compared to earlier cycles (cycles 12-17), and the differ- 
ences in their means are statistically important. Thus, cycles of late have been more robust than 
earlier cycles and the higher parametric values may result from a long-term secular increase over 
time. 16 ’ 17 

Figure 5 shows the cyclic variation of specific timing parameters, where 7(1) is the elapsed 
time in months from E(Rm ) to the epoch of ARmaxpos (£lA/bnaxpos)), t( 2) is the elapsed time 
in months from E(ARmaxpos) to E(RM), t( 3) is the elapsed time in months from E(RM) to the 
epoch of T/huaxneg (£fA/?maxneg)), and t( 4) is the elapsed time in months from E(ARmaxneg) to 
E(Rm) of the next cycle. Figure 6 shows the same timing parameters but using the 12-mma values of 
T/huaxpos and ARmaxneg. In both figures, the horizontal lines are the means (given to the right in 
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Figure 2. Variation of AR for elapsed time in months / = 0-132 mo from E(Rm). 
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Figure 3. Variation of 12-mma of AR for elapsed time in months t- 0-126 mo from E(Rm). 
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Figure 4. Variation of parametric values over cycles 12-23 


Sunspot Cycle Number Sunspot Cycle Number Sunspot Cycle Number 


Based on AR 




(b) 



mean = 25.8 
sd = 9.3 


mean = 18.8 
sd= 13.5 



Figure 5. Variation of specific timing signatures based on AR. 
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Based on 12-mma of AR 


(a) 



mean = 51.4 
sd = 9.6 



mean = 30.2 
sd = 6.8 


(C) 



mean = 23.0 
sd= 11.4 



Figure 6. Variation of specific timing signatures based on 12-mma of AR. 
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each subpanel along with the sd). No significant differences in the timing parameters are apparent, 
comparing earlier cycles with those of late. Thus, on average, the greatest positive rate of growth in 
R occurs about 2 yr after Rm and about 1.5-2 yr before RM, and the greatest negative rate of growth 
in R occurs about 2-2.5 yr after RM and about 4-5 yr before succeeding cycle Rm. Table 2 provides 
a convenient summary of the cyclic parametric values and their occurrences relative to E(Rm). 


Table 2. Selected parametric values for sunspot cycles 12-23. 


Cycle 

Rm(E(Rm)) 

RM(E(RM)) 

ASC 

PER 

4Rmaxpos(f) 

aRmaxneg(f) 

12-mmaof 

/iRmaxpos(t) 

12-mmaof 

4Rmaxneg(f) 

12 

2.2 (12-1878) 

74.6(12-1883) 

60 

135 

4.1 (14) 

-4.1 (83) 

2.47(19) 

-2.55 (88) 

13 

5.0 (03-1890) 

87.9(01-1894) 

46 

142 

4.6 (16) 

-3.6 (57) 

3.36(14) 

-1.70 (69) 

14 

2.6 (01-1902) 

64.2 (02-1906) 

49 

139 

3.9 (39) 

-5.4 (51) 

2.00 (20) 

-2.08 (96) 

15 

1.5 (08-1913) 

105.4 (08-1917) 

48 

120 

8.5 (42) 

-5.3 (75) 

3.82(41) 

-2.69 (72) 

16 

5.6 (08-1923) 

78.1 (04-1928) 

56 

121 

6.3 (22) 

-5.5 (82) 

3.07 (23) 

-2.70 (78) 

17 

3.4 (09-1933) 

119.2 (04-1937) 

43 

125 

6.4 (39) 

-4.7 (73) 

4.25 (35) 

-2.47 (73) 

18 

7.7 (02-1944) 

151.8(05-1947) 

39 

122 

8.6 (33) 

-6.9 (73) 

5.95 (31) 

-4.30 (71) 

19 

3.4 (04-1957) 

201.3(03-1958) 

47 

126 

10.8 (22) 

-6.2 (75) 

8.69 (21) 

-4.79 (79) 

20 

9.6 (10-1964) 

110.6 (11-1968) 

49 

140 

6.4 (22) 

-5.3 (73) 

3.97(19) 

-3.21 (74) 

21 

12.2 (06-1976) 

164.5(12-1979) 

42 

123 

8.1 (24) 

-8.4 (74) 

5.83 (20) 

-4.02 (78) 

22 

12.3 (09-1986) 

158.5(07-1989) 

34 

116 

10.6 (21) 

-8.3 (64) 

6.92 (21) 

-5.45 (66) 

23 

8.0 (05-1996) 

120.8 (04-2000) 

47 

- 

5.5 (40) 

-5.3 (77) 

3.57 (20) 

-3.43 (78) 


Recall from figure 4 that there is the hint that cycles of late have been more robust than 
earlier cycles and that this behavior may be the result of a long-term secular increase over time. 
Table 3 provides a convenient summary of the statistics, comparing parametric values against sun- 
spot cycle number. Indeed, all parameters except PER appear to correlate well against sunspot cycle 
number. Table 3 gives the coefficient of linear correlation (r), the coefficient of determination (r 2 ) 
(a measure of the amount of variance explained by the regression), the v-intercept in the regression 
equation (a), the slope of the inferred regression (b), the standard error of estimate (se), and the con- 
fidence level (cl), where cl >95 percent is considered statistically important (cl >90 percent is of mar- 
ginal statistical importance). Projections (90-percent prediction intervals) for cycle 24 based on the 
inferred regressions are also given. Hence, before cycle 24 has officially started, presuming the validity 
of the inferred regressions, cycle 24 is expected to have Rm = 1 1 .5 ± 4.5 (this will be exceeded low, for 
R measured 5.9 in September 2007), RM= 166.9±64.1, ASC= 39.3 ± 10.7 mo, ARmaxpos = 9.4±3.7, 
ARmaxneg = -7.6±2.0, 12-mmaARmaxpos = 6.55 ±3.03, and 12-mmaARmaxneg = -4.83± 1.43. Pre- 
suming cycle 23 has a PER = 141 mo (implying E(Rm) for cycle 24 in March 2008), PER for cycle 24 
is estimated to be about 125 ±18 mo. 


2.2 Correlative Results 


2.2.1 ARmaxpos 

The role of ARmaxpos as a predictor for RM will be examined in this section. Figure 7 shows 
the scatter plot of RM versus ARmaxpos, where TRmaxpos is recalled as the monthly TRmaxpos 
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Table 3. Correlations against sunspot cycle number. 


Parameter 

r 

r 2 

a 

b 

se 

cl (%) 

Cycle 24 
(Predicted)** 

Rm 

0.780 

0.608 

-8.248 

0.821 

2.504 

>99.5 

11.5±4.5 

RM 

0.594 

0.353 

20.868 

6.085 

35.375 

>95 

1 66.9 ± 64. 1 

ASC 

-0.591 

0.349 

66.614 

-1.140 

5.902 

>95 

39. 3± 10.7 

PER * 

-0.510 

0.260 

152.045 

-1.409 

8.306 

<90 

118.2 ± 15.2 

ztRmaxpos 

0.558 

0.312 

0.620 

0.364 

2.029 

>90 

9.4±3.7 

ztRmaxneg 

-0.685 

0.470 

-0.806 

-0.283 

1.108 

>98 

-7.6 + 2. 0 

12-mmaof aRmaxpos 

0.578 

0.334 

-1.031 

0.316 

1.674 

>95 

6.55±3.03 

12-mmaof aRmaxneg 

-0.746 

0.556 

0.861 

-0.237 

0.789 

>99.5 

-4.83 ±1.43 


’Excludes cycle 23 (presuming PER= 141 for cycle 23, then the numbers are, respectively, -0.223, 0.050, 139.446, -0.587, 9.663, <90, and 125.3 ± 17.5) 
**90-percent prediction interval 


(growth in R) during the ascending portion of the sunspot cycle (the ascending inflection point). The 
number beside each of the filled circles identifies the specific sunspot cycle. As an example, cycle 19 
(the uppermost filled circle) had its greatest positive change in R (= 10.8) between t- 22 (R = 98.5) 
and t = 23 (7? = 109.3) and it is this number that has been correlated with its later occurring RM 
(= 201.3). The thick straight line (denoted y) is the inferred regression line and the thin vertical and 
horizontal lines are the medians for the two parameters. The inferred regression is computed as 
y = 11.365 + 15. 519.x; and has r = 0.867, yielding r 2 = 0.755 (meaning that about three-fourths of the 
variance in RM can be explained by the inferred regression based on the behavior of ARi naxpos), 
and se = 21.8. The regression is found to be statistically significant at better than the 0. 1-percent level 
of significance (or cl >99.9 percent). The result of Fisher’s exact test for 2x2 contingency tables is 
also given, 18 which indicates that the probability (P) of obtaining the observed result, or one more 
suggestive of a departure from independence (chance), is P =12.1 percent. It is apparent then that 
by monitoring the month-to-month rate of change in R, the later occurring RM can be continually 
estimated. For ARmaxpos <6.4, it is anticipated that RM <121; for ARmaxpos >6.4, it is anticipated 
that RM >105. 

Instead of using T/dnaxpos as the estimator for the later occurring RM. RM can be compared 
against the slope at Tidnaxpos, where the slope is computed as R at ZfiA/bnaxpos) minus Rrn divided 
by the elapsed time in months between E{AR\ naxpos) and E(Rm). Figure 8 displays the scatter plot 
of RM versus the slope at FfARmaxpos). As before, the number beside the filled circles identifies the 
sunspot cycle number, the thick line is the inferred regression, the thin lines are the medians, and the 
results of linear regression analysis and Fisher’s exact test for the 2x2 contingency table are given. As 
an example, cycle 19 (the uppermost filled circle) had ARmaxpos= 10.8 occurring at t = 22 mo. The 
slope for cycle 19 is then computed to be (98.5 - 3.4)/22 = 4.32, and it is this number that is correlated 
with the later-occurring RM (= 201.3). 

Clearly, figure 8 provides a much-improved prediction of RM as compared to using fig- 
ure 7 because of the large reduction in se (reduced by half, from 21.8 to 10.6), but not knowing 
exactly when ARi naxpos occurs presents a problem. For example, cycles 14, 15, 17, and 23 had 
their actual Zifhnaxpos values very late in their ascents (t = 39, 42, 39, and 40 mo, respectively), 
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Figure 7. Scatter plot of RM versus Ztifrnaxpos. 
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Figure 8. Scatter plot of RM versus slope at FfARmaxpos). 
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although smaller local peaks occurred earlier in their ascents (3.2 at t - 16 and 3.5 at t = 30 for cycle 
14, 5.4 at t- 17 for cycle 15, 4.8 at t = 24 and 5.4 at t- 33 for cycle 17, and 5.2 at i = 20 for cycle 23). 
Using these earlier values, cycle 14’s RM would have been predicted to be about 61 (at t- 16) and 
66 (at t = 30), both values very close to its actual RM= 64; cycle 15’s RM would have been predicted 
to be about 95 (at t- 17), close to but below its actual RM- 105; cycle 17’s RM would have been 
predicted to be about 86 (at t = 24) and 95 (at ^ = 33), both close to but below its actual RM= 119; 
and cycle 23’s RM would have been predicted to be about 92 (at I - 20), close to but below its actual 
RM- 121. Because the 90-percent prediction intervals associated with each predicted RM is ±39.5, 
all of the actual RM values would have fallen within the prediction interval. 

Instead, using the estimated slope at the time of the earlier peaks, it would have been pre- 
dicted that cycle 14’s RM to be about 60 (at t- 16) and 72 (at ? = 30), cycle 15’s RM to be about 
84 (at t- 17), cycle 17’s RM to be about 91 (at t- 24), and cycle 23’s RM to be about 91 (at t- 20). 
Because the 90-percent prediction intervals associated with each predicted RM is± 19.3, only the 
actual RM values for cycles 14 and 15 would have fallen within the prediction interval; cycles 17 
and 23 would have fallen just outside the upper limit of the prediction interval (by about 10 units of 
sunspot number). 

Because many of the cycles either had their Tifrnaxpos or alternative early cycle local peaks 
before t = 30, it may be that RM can be predicted on the basis of the collective behavior of Tihnaxpos 
in comparison to RM before ? = 30 mo. Figure 9 shows the variation of the r 2 s determined month- 
by-month from E(Rm) based on a comparison of RM and the month-to-month change in R(AR)(t) 
for t = 0-26 mo. It is found that the peak r 2 occurs at t = 25 mo; although, as early as t = 14 mo, AR(t) 
provides a statistically meaningful estimate for RM (the earliest estimate appears to be at t- 10 mo, 
but it has a rather large se, equal to ± 34 mo). In figure 9, the cross-hatched pattern signifies those /s 
when a statistically meaningful correlation is inferred between RM and AR(t). 

Figure 10 displays the scatter plots of RM versus AR(t= 10) (panel (a)) and AR(t = 25) (panel 
(b)). Clearly, by monitoring AR(t), an increasingly accurate estimate for RM can be effected, espe- 
cially at t = 25 mo. Such an approach would have suggested its RM would likely be about 105.6 ± 33.6 
(the 90-percent prediction interval) for cycle 23. 

2.2.2 12-mma of AR 

The role of the 12-mma of AR as a predictor for RM will be examined in this section. 
Figure 11 shows the scatter plot of RM versus the maximum positive value of the 12-mma of AR. 
The structure of the chart follows that used in figures 7, 8, and 10. Figure 11 is strikingly similar 
to figure 8, both having r = 0.97, r 2 = 0.94, and se- 10.6. Thus, by monitoring the 12-mma of AR, 
increasingly accurate predictions of the later occurring RM can be effected. The disadvantage of the 
12-mma of AR(t) as compared to AR(t) is that it lags that of AR(t) by several months, so its use as a 
predictor of RM appears more appropriately to be that of a confirmatory role. 

Figures 12 and 13 are equivalents to figures 9 and 10, but based on the comparison of RM 
versus the 12-mma of AR(t). From figure 12, it is found that as early as ? = 8 mo, the value of the 
12-mma of AR(t) provides a statistically meaningful estimate for RM, with the best predictor 
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t, Elapsed time in months from E(Rm) 


Figure 9. Variation of r 2 (AR(t)). 


occurring at t = 22 mo. For cycle 23, such an approach would have suggested its RM likely would be 

1 1 1.3 ± 25.8 (the 90-percent prediction interval). 

2.2.3 A/?maxneg 

The role of ARmaxneg as a predictor of PER for the current cycle and of Rrn and RM for 
the following cycle will be examined in this section. Figure 14 displays the scatter plot of PER versus 
ARmaxneg. The structure of the chart follows that of previous charts. Based on Fisher’s exact test for 
2x2 contingency tables, it is found that the P of obtaining the observed result, or one more suggestive 
of a departure from independence (chance), is P= 17.5 percent. Based on linear regression analysis, 
a statistically important correlation is inferred to exist between PER and ARmaxneg, having r = 0.66, 
r 2 = 0.44, se-1 .2, and cl >95 percent. From the known value of ARmaxneg for cycle 23 (= -5.3 at 
t-11 mo or 30 mo past E(RM)), it is inferred that PER for cycle 23 should be about 130 ±13 mo 
(the 90-percent prediction interval). Thus, on the basis of the validity of the inferred regression, it is 
noted that there is only a 5-percent chance that cycle 23 will persist longer than 143 mo, indicating 
E(Rm) for cycle 24 should be expected before May 2008. 
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Figure 10. Scatter plots of RM versus AR{t = 10) and RM versus AR(t = 25). 
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Figure 11. Scatter plot of RM versus (12-mma of AR)maxpos. 



(, Elapsed time in months from E(Rm) 

Figure 12. Variation of r 2 (12-mma of AR(t)). 

Figures 15 and 16 show scatter plots of Rm(n + 1) and RM(n + 1) versus ARmaxneg, respec- 
tively. Neither plot is inferred to be statistically important. Thus, while T/huaxneg seems to provide 
an indication for the expected length of an ongoing cycle, it does not provide any indication as to the 
expected size of the Rms and RMs for the following cycle. 

2.2.4 (12-mma of Aft)maxneg 

The role of (12-mma of Ti?)maxneg as a predictor of PER for the current cycle and of Rm 
and RM for the following cycle will be examined in this section. Figure 17 displays the scatter plot 
of PER versus (12-mma of AR)maxneg. As found for PER versus ARmaxneg (fig. 14), a statisti- 
cally important correlation is inferred between PER and (12-mma of AR) maxneg, having r = 0.63, 
r 2 = 0.39, se = 1.5, and cl >95 percent. From the known value of (12-mma of T/?)maxneg for cycle 23 
(= -3.43 at / = 78 mo or 31 mo past E(RM)), it is inferred that PER for cycle 23 will be about 127 ± 14 
mo (the 90-percent prediction interval). Thus, on the basis of the inferred regression, there is only a 
5-percent chance that cycle 23 will persist longer than 141 mo, indicating E(Rm) for cycle 24 should 
be expected before April 2008. 


19 


o 



Figure 13. Scatter plots of RM versus 12-mma of AR(t = 8) and RM 
versus 12-mma of Z\R(t = 22). 
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Figure 14. Scatter plot of PER versus Tihnaxneg. 


Figures 18 and 19 show scatter plots of Rm(n+ 1) and RM(n + 1) versus (12-mma of AR) 
maxneg, respectively. Neither plot is inferred to be statistically important. Thus, while (12-mma of 
AR)maxneg seems to provide an indication for the expected length of an ongoing cycle, it does not 
provide any indication as to the size of the Rms and RMs of the following cycle. 

2.2.5 AR(T) and 12-mma of AR(T) 

The roles of AR(T) and 12-mma of AR(T) for elapsed times in months from E(RM), T, as 
predictors of PER, Rm(n + 1 ), and RM(n + 1) will be examined in this section. Figure 20 shows 
the scatter plots of PER versus AR(T=30) (panel (a)) and PER versus 12-mma of AR(T= 31) 
(panel (b)). The two plots represent the best fits of PER versus AR(T) and PER versus 12-mma 
of AR(T) for T-0 40 mo. Apparently, the reason for the success of these inferred regressions is 
the fact that AR(T) and the 12-mma of AR(T) are near their maximum negative values at these 
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Figure 15. Scatter plot of Rm ( cycle n + 1) versus Zlihnaxneg. 


elapsed times relative to E(RM). Both fits are statistically important. Based on the known value 
of AR(T= 30) during the descending portion of cycle 23 (= -5.3), it is inferred that its PER will be 
about 122 ± 14 mo (the 90-percent prediction interval). Based on the known value of the 12-mma 
of AR(T= 31) during the descending portion of cycle 23 (= -3.43), it is inferred that its PER will 
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Figure 16. Scatter plot of RM(cyc\e n + 1) versus ARmaxneg. 
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Figure 17. Scatter plot of PER versus (12-mma of T/?)maxneg. 


be about 125 ± 12 mo (the 90-percent prediction interval). Since cycle 23 has already persisted for 
135 mo, E(Rm) should be expected very soon. However, because new cycle spots typically occur 
either simultaneously with up to several months prior to E(Rm) and because the first confirmed new 
cycle spot for cycle 24 occurred in January 2008, it appears that E{Rm) for cycle 24 will occur out- 
side these 90-percent prediction intervals. (The interval from May 1996- August 2007 corresponds to 
135 mo. Because cycle 23 undoubtedly will have a PER more like cycles 12-14 and 20 rather than 
cycles 15-19, 21, and 22, inclusion of cycle 23’s PER will greatly weaken the inferred regressions. 
Hence, the inferred correlations are probably a fluke, with cycles more likely distributed preferen- 
tially as short- and long-period cycles, rather than following the inferred linear regression lines.) 

Figure 21 shows the scatter plots of Rm(n+ 1) versus AR(T =31) (panel (a)) and Rm(n + 1) 
versus 12-mma of AR(T=21) (panel (b)), the most statistically important correlations. While 
the first plot (panel (a)) is only of marginal statistical importance (cl >90 percent), the latter one 
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Figure 19. Scatter plot of RM(cycle n + 1) versus (12-mma of Z\i?)maxneg. 
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Figure 20. Scatter plots of PER versus AR(T = 30) and 12-mma of AR(T = 31). 
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Figure 21. Scatter plots of Rm { cycle n + 1) versus AR(T = 31) and 12-mma of AR(T = 27). 
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(panel (b)) is found to be highly statistically important (cl >99.5 percent). Based on the known value 
(= -3.2) of cycle 23’s AR(T =31), the Rm for cycle 24 is estimated to be about 7 ± 6 (the 90-percent 
prediction interval); while, based on the known value (= —2.88) of cycle 23’s 12-mma of AR(T =27), 
the Rm for cycle 24 is estimated to be about 8 ±4 (the 90-percent prediction interval). Based on the 
2x2 contingency table, Rm(24) can be expected to be >5.6. The lowest value of the 12-mma of R 
observed to date during the declining portion of cycle 23 measures 5.9 (September 2007), so obvi- 
ously E(Rm ) for cycle 24 is most imminent. 

Attempts to find a statistically important correlation between RM(n + 1) and AR(T) and 
between RM(n+ 1) and the 12-mma of AR(T) proved fruitless. However, for completeness sake, fig- 
ure 22 is included, which shows the scatter plots of RM(n+ 1) versus AR(T= 14) (panel (a)) and 
RM(n+ 1) versus the 12-mma of AR(T=44) (panel (b)). Based on the known value (= 1.9) of cycle 
23’s AR(T= 14), the RM for cycle 24 is estimated to be about 139 ±67 (the 90-percent prediction 
interval); while, based on the known value (= -1.88) of cycle 23’s 12-mma of AR(T=44), the RM of 
cycle 24 is estimated to be about 137 ± 69). 




Figure 22. Scatter plots of RM( cycle n + 1) versus AR(T = 14) and 12-mma of AR(T = 44). 


29 


2.2.6 Slope ASC and Slope^g 

Figure 23 displays the scatter plots of Slope^^ (panel (c)) and Slope^,^ (panel (a)) versus 
ARmaxpos and Slope /1<iC (panel (d)) and Slope^,^ (panel (b)) versus (12-mma of AR)maxpos, where 
Slope^^ is dehned as (RM Rm)/ASC and Slope^,^ is defined as ((Rm for cycle n + 1 )-(RM for 
cycle n))/DES for cycle n. All correlations are inferred to be highly statistically important, especially 
Slope^,^ versus (12-mma of AR)maxpos. Estimates (90-percent prediction intervals) of Slope^,^ 
for cycle 23 are made based on the observed values of ARmaxpos and (12-mma of AR)maxpos. 

Figure 24 shows scatter plots of Slope^,^ versus ARi naxneg (panel (a)), (12-mma of AR) 
maxneg (panel (b)), and Slope^^ (panel (c)). Again, all correlations are inferred to be highly statisti- 
cally important, especially Slope DHS versus Slope /J ()C . Estimates (90-percent prediction intervals) of 
Slope Df ,- S for cycle 23 are made based on the observed values of ARmaxneg, (12-mma of AR)maxneg, 
and Slope^ AC . 


Of the various estimates for SI ope the best estimate (smallest 90-percent prediction inter- 
val) appears to be the one based on (12-mma of AR)maxpos. Hence, cycle 23’s Slope /)/;A likely will 
measure about -1 . 1 5 ± 0.29. 

2.2.7 Cycle 23 Behavior: May 1996-August 2007 

Figure 25 displays the 12-mma of R (panel (b)) for elapsed times t = 0 (May 1 996)— 1 35 (August 
2007). The 12-mma of AR for elapsed times t- 0 (May 1996)— 128 (January 2007) is also shown. 
For convenience, the values of Rm, RM, ASC, DES, PER, (12-mma of AR)maxpos, and (12-mma 
of AR)maxneg and the dates of occurrence of E(Rm), E(RM), £(( 12-mma of A7?)maxpos), and 
£((12-mma of Ti?)maxneg) are given. Thus, RM for cycle 23 measured 120.8 (occurring at / = 47, 
April 2000) and the value of R at /= 135 mo (August 2007) measured 6.1. The current Slope^,^ is 
(6.1 - 120. 8)/88 =-1.30, which is within the window of the estimated value for cycle 23’s Slope^,^ 
(= -1.15 ±0.29), inferring that E(Rm ) for cycle 24 is imminent. Recall from figure 6 that the time 
(l(4)) from £((12-mma of TR)maxneg) to E(Rm) of the following cycle averages about 51 mo (with 
sd= 10), inferring that E(Rm ) for cycle 24 should follow November 2002 (l = 78) by about 51 ± 18 mo 
(the 90-percent prediction interval), or that there is a 95-percent P that cycle 24 will have its official 
onset before September 2008. (Since R is now known through September 2007, the actual Slope^,^ 
is -1.29.) 
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Figure 23. Scatter plots of Slope^^ versus ARmaxpos and (12-mma of AR)maxpos 
and of Slope^,^ versus Z\i?maxpos and (12-mma of AR)maxpos. 


31 



' o 

CO 

*=* 

CD 

Q. 

O 

CO 



Figure 24. Scatter plots of Slope^,^ versus z\i?maxpos, (12-mma of AR)maxpos, and Slope^ 5C 
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tion of cycle 23’s R and 12-mma of AR for t = 0-135 (May 1996-August 2007). 
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3. SUMMARY 


Various methods, 4 ’ 19-22 those using geomagnetic precursor information in particular, have 
been used to successfully predict RM several years in advance of its occurrence. For cycle 24, the 
next sunspot cycle, precursor techniques suggest that it will have an RM equal to about 1 30 ±14 (see 
Wilson and Hathaway 22 and references contained therein), peaking about 44 ± 5 mo after E(Rm), or 
about November 201 1 (± 5 mo) if the official start of cycle 24 is, indeed, March 2008, as predicted 
by the NOAA Solar Cycle 24 Prediction Panel. 14 This prediction for the size of cycle 24’s RM is 
more supportive of the higher consensus prediction of the NOAA panel (140 ±20) than its lower 
consensus prediction (90 ± 10), where the two consensus predictions arise from consequences of two 
different dynamo-related techniques. So, it is with great anticipation 23 that solar researchers await 
the start and rise of cycle 24. 

Once a cycle officially gets underway, predictions of the level of future activity based on 
autoregressive techniques 24 ’ 25 and other curve-fitting methods 6 ’ 9 ’ 19 can be employed. As found in 
this Technical Publication and in previous studies, 2 ’ 6 ’ 7 ’ 9 ’ 19 the size ( RM) of the cycle can be better 
predicted about 2 yr into the rise of a sunspot cycle, especially once its ascending inflection point has 
been clearly discerned. 

In this Technical Publication, it has been noted that, on the basis of cycles 12-23 (the most 
reliably determined sunspot cycles), cycles can be loosely grouped into three groups based on RM 
or two groups based on cycle duration. Large-amplitude cycles (cycles 18, 19, 21, and 22) tend to 
have higher values of Rm, rise to RM more quickly (shorter ASQ, and have shorter cycle duration 
(shorter PER) than either average-amplitude cycles (cycles 15, 17, 20, and 23) or small-amplitude 
cycles (cycles 12, 13, 14, and 16). Likewise, short-period cycles (cycles 15-19, 21, and 22) tend to rise 
more quickly to RM and to have higher values of Rm and RM than long-period cycles (cycles 12-14, 
20, and 23). It is also noted that, over the course of cycles 12-23, there have been statistically impor- 
tant secular rises in Rm and RM and in the values of R at the ascending (more positive) and descend- 
ing (more negative) inflection points. On the basis of the secular increases, estimates have been made 
for cycle 24; namely, its Rm should measure about 1 1 .5 ± 4.5, its RM about 167 ± 64, its ASC about 
39 ±11 mo, the value of its ARmaxpos about 9.4 ±3.7, the value of its ARmaxneg about -7.6 ±2.0, 
the value of its 12-mmaARmaxpos about 6.55 ±3.03, and the value of its 12-mmaARmaxneg about 
-4.83 ± 1.43. Its PER will measure about 125 ±18 mo. 

This study has shown that the later occurring RM can be continually estimated by moni- 
toring the month-to-month rate of change in R. For values of ARmaxpos <6.4, it is expected that 
RM <121, while for values of ARmaxpos >6.4 it is expected that RM >105. Based on the slope deter- 
mined at ARmaxpos, RM can be predicted with even higher precision (se = 10.6). It was found that, 
as early as 10 mo after E(Rm), the value of AR(t) can be used to estimate the later occurring RM, 
although the best result is determined about 2 yr into the cycle. Use of the 12-mma of AR to predict 
RM provides similar accuracy as compared to using the slope at ARmaxpos. Values of the (12-mma 
of AR )maxpos <4.25 strongly suggest RM <121, while higher values suggest a larger RM. 
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Both ARmaxneg and its 12-mma value appear to correlate with PER, such that cycle 23’s 
PER is expected to persist no longer than 143 mo, inferring that cycle 24 should have its official start 
no later than May 2008. An examination of AR and its 12-mma value relative to E(RM) suggests that 
Rm for cycle 24 will measure about 7.6 ± 4.4. R measured 5.9 in September 2007, so Rm for cycle 24 is 
expected very soon. The (12-mma of AR)maxpos suggests that the Slope^,^ for cycle 23 will measure 
about -1.15 ±0.29. The computed current Slope^,^ for cycle 23 is -1.29 through September 2007 
and is shrinking with the passage of time. On average, the length of time in months from (12-mma of 
AR)maxneg to Rm is 51 ± 18 mo (the 90-percent prediction interval). Thus, there is only a 5-percent 
chance that E(Rm) for cycle 24 will occur after August 2008. 
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